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DESCRIPTION 

METHOD OF POSmON DETERMINATION 

The present invention relates to a method of position determination 
in a radio system and to apparatus for position determination 1 n a radio 
system. 

European Patent Publication 0 654 677 A describes a method of 
processing a signal in a multi-path environment. 

However, such methods may not perfomri particularly well In the 
presence of diffuse reflections especially when such components are very 
plentiful and closely spaced, as for example occur in indoor environments 
involving rough uneven surfaces and clutter. Accordingly, Imprecise 
estimations and hence Inaccurate positioning result. 

An object of the present invention Is the avoidance or elimination of 
diffuse multi-path inaccuracies. 

Another object of the present invention is to provide good 
positioning accuracy. 

Another object of the present Invention Is to remove diffuse 
scattering components from received signals. 

According to the present Invention, there i s p rovided a method of 
position determination In a radio system, the method comprising 
correlating a signal received at a unit with a replica signal at the unit, and 
processing the correlated signal with an optimisation function comprising 
an exponential term In combination with a second term. 

In tills way, the method of the present Invention may provide 
accurate position measurement In indoor environments. 

Preferably, the exponential term is in the form Be"°*, and ihe second 
term is of the form: 
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The method^ may comprise effecting an integration witli tine replica 
signal, and/or fitting the optimisation function and a Line-of Sight 
correlation function with a set of parameters, and/or superposing the 
diffuse correlation output with a Line-of-Sight function output and fitting 
5 with correlation data of known values for the Line-of-Sight output. 

The method may comprise first operating a multipath mitigation 
technique to effect correlation of the received and replica signals. 

The multipath mitigation technique may comprise a Multipath 
Estimating Delay Locks Loop (IVIEDLL) technique, or a Minimum Mean 
•.<ib Square Error (MMSE) technique. 

According to the present invention, there is also provided a 
computer program product directly loadable into the internal memory of a 
digital computer, comprising software code portions for performing the 
method of the present invention when said product is run on a computer. 
15 According to the present invention, there is also provided a 

computer program directly loadable i nto t he i nternal memory of a digital 
computer, comprising software code portions for performing the method of 
the present invention when said program is run on a computer. 

According to the present invention, there is also provided a camer, 
20 which may comprise electronic signals, for a computer program 
embodying the present invention, and/or electronic distribution of a 
computer program product, or a computer program, or a carrier of the 
present invention. 

According to the present invention, there is also provided apparatus 
25 for position determination of a radio system, the apparatus comprising 
means to correlate a signal received at a unit with a replica signal at the 
unit, and means to process the correlated signal with an optimisation 
function comprising an exponential term in combination with a second 
term. 

30 The present invention is applicable to radio systems especially 

involving GPS and GSM networks, particulariy in indoor environments and 
can be implemented in hardware and/or software. 
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In order that the present invention may more readily be understood, 
a description is now given, by way of example only, reference being made 
to the accompanying drawings in which:- 

Figure 1 is a schematic d iagram of a receiver unit of the present 
invention; 

Figure 2 is a hardware implementation of the unit of Figure 1; 
Figure 3 is a flow diagram of the signal processing operation of the 
method of the present invention; and 

Figure 4 depicts the surfaces that diffusers lie on. 

Figure 1 shows a mobile radio receiver unit 1 of the present 
invention which is located in a room 2 with a variety of objects including 
chairs 3, desks 4, bookcase 5 and reflective surfaces e.g. walls 6. These 
objects act in combination to generate a large number of diffuse reflections 
for GPS signals of the type for reception at unit 1 . 

Details of the componeints in a hardware Implementation of unit 1 
are shown in Figure 2, having an antenna 10 to receive incoming signal 
R(t) and pass It to multiplier block 1 1 for multiplication with output C from 
the canier generator 12 in order to mix down signal R(t). The resultant 
signal Is passed to a block of correlators 1 3 for sequential correlation of 
the mixed down signal and the replica code from code generator, in order 
to produce correlation data outputs from the I (Imaginary) and Q (real) 
channels using Input signals from code generator 14. 

The correlation data is then passed to signal processing block 1 5 
which adjusts appropriately both the code generator 13 and carrier 
generator 12. 

Figure 3 Is a flow diagram of the major steps in the correlation 
operation for signal processing to remove the diffuse components as in the 
present invention. 

Step 1 of the flow diagram involves standard signal processing for 
position measurement In the presence of multi-path signals, typically 
accomplished by the Multipath Estimating Delay Lock Loop (MEDLL) 
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method generally as described in European Patent Publication No. EP 0 
654 677 A. 

Step 1 involves using this method to provide an optimal fit using N 
"triangles" thereby yielding estimates for the multipath parameters:- 

{a_0.....a_N} 
{tJ,....t^N} 
{q_1.....t_N} 

Then step 2 adds to the standard conrelation operation by using the 
initial value for t_0 obtained from the MEDLL method, incorporating small 
initial values for B and a, and then correlating the function 

f(r) = 5a-(r,^/l-^ + l) (1) 

with the replica signal s(t-r) . 

The correlation is performed using multiplier block 1 1 so that the 
received signal r(t) is multiplied by the replica code on the carrier wave to 
mix down the received signal as well as obtain the correlation data. This is 
then summed by integrator 1 6 in Step 3 to give correlation function data for 
the I and Q channels, the congelation function being complex due to the 
phase of the signal. 

In step 4, f (r) is fit with the LIne-of-Sight triangle with parameters 

{a_0, t_0, clP] 
to the correlation data using l<nown values 
{a_1....a_N} 
{t„1,....t_N} 
{qJ....CLN} 

where 

ao,...aN are the amplitudes 
ro,....T^ are the delays 

qo,....qN are the phases 
of all the reflected components in the received signals. 
Step 5 involves finding values for 

B. 
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a, 

a_0, 
clO. 
LO 

and E (the mean square en^or), 

trying different values until a global minimum error Is reached, and 
storing all the derived values. 

Step 6 involves integrating the same operations for the respective 
values, and then feeding them back to Step 2 to reiterate the process. 

Thus the function f(r) which Is characterised in equation (1) acts 
efficiently to clean the correlation functions of diffuse reflections. 

As a variant, a software implementation of the correlation operation 
incorporates an algorithm to effect: 

• performing the correlation of the local replica code with the received 
signal data using multiplication and summation of the sampled values 
of the local code and the received signal^^ to produce demodulated 
correlation data; 

• operating the conventional p osition m easurement method by fitting N 
amount of triangles of fixed width 2* chip width to the demodulated 
correlation data such thatN triangles give the minimum g lobal error; 
and 

• operating the f (r) of equation 1 and Stage 2. 

In greater detail, Figure 3 represents an additional stage of signal 
processing, in accordance with the present invention, for the purpose of 
cleaning a received correlation function of remaining diffuse reflection 
components for a number of multipath mitigation techniques, for example a 
Muttipath Estimating Delay Lock Loop (MEDLL) or a Minimum Mean 
Square Error (MMSE) technique. 

Although the method shown in Figure 3 and the following 
description both refer to a f^EDLL technique, it should be understood that 
the improvement to this technique described hereinafter is equally 
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applicable to other forms of multipath mitigation technique, one example of 
which is a technique referred to as the Minimum Mean Square Error 
(MMSE) method. 

The additional stage in the signal processing utilises a cleaning 
function, /(r), in the maximum likelihood estimator (MEDLL) or Bayesian 
estimator (MMSE), making possible the removal of the diffuse scattered 
components from the received correlation function data. 

The present invention allows for the optimisation of the 
implementation of a maximum likelihood estimator, using a novel 
correlation function that is a sum. of a line-of-sight (LOS) correlation 
function and a diffuse correlation function, both of which are dependent on 
the LOS delay Tq . 

The improved estimator utilises an approximation of the diffuse 
scattering process. 

The maximum likelihood estimator of a conventional MEDLL 
technique as currently in use is given by: 

L ^ lj[D(T)-f^a„e''"P{T„)fxlT * (3) 

V «*=o 

where 

R(r)^f^a„e^'-PiT„) . (4) 

and where P(t„) is a triangle of unit height centred at time r = r„ 
with a base of two chip widths 21^ , 

phase of the nth reflection , 
nth amplitOde a„ and 
nth delay r„ 

which are parameters in the correlation function model. 
The term D(t) , where: 

D(T)=jrit)s(t-'T)jdt , (5) 
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results from the received signal data, r(t), correlated with the 
replica signal, sit) , in the receiver to obtain the con-elation function data. 

Thus, the maximum likelihood estimator for the multipath 
parameters {aQ..,aj^,0Q..j0j^,rQ,.s2^} used in the MEDLL m ethod i s g iven 

by those parameter values that minimise the estimator L . 

In deriving the f(r) function of equation (1), three main factors were 
considered, namely: 

i) the combined amplitude of randomly-scattered superposed rays as 
path length or time increases; 

ii) the number of randomiy-piaced diffusers as path iength or time 
increases; and 

iii) " the free-space decay of ray amplitudes as path length or time 

increases. 

In order to derive the form for the randomly-scattered superposed 
rays, mentioned in i) above, the diffuse reflectors are assumed to be 
homogeneously distributed since they can be viewed as scattering the 
rays from randomly placed locations. Moreover, at any given path iength, 
all ray amplitudes are approximately equal, and so the random scattering 
is only manifested in the phases of individual received rays. 

As time increases, the likelihood of any two received rays being in 
phase with one another is reduced, as rays with longer reflected paths are 
scattered further. Hence, the phases of individual rays of the same path 
length become increasingly divergent, and consequently, the amplitudes of 
the individual rays are more likely to cancel one another out upon 
superposition. 

The exponential, decay of the random-scattering effect is governed 
by the equation: 

f(t) = Ae-''' . (6) 
where / is the propagation delay of the rays. The constant b 
determines how quickly the .function tends to zero. Since more re-radiatlon 
occurs at shorter wavelengths, the wavelength, X , has the property: 
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Replacing b with tiie physical parameter X yields the following form 
for the random phase approximation of the received signal: 

/(0 = ^e-*'^^ (7) 

where A is the amplitude of the transmitted signal. 

In order to determine the number of diffusers with respect to time, 
mentioned in ii) above, it is assumed that each diffuser is a point 
homogeneously distributed in space. For any fixed path length for a signal 
transmitted from a transmitter, y4, to a receiver, B , the reflectors must lie 
on the surface of an ellipsoid. Therefore, t he number of reflectors N is 
proportional to the surface area of an ellipsoid with foci separated by the 
distance AB . 

Figure 4 shows the surfaces that diffusers lie on for a fixed path 
length r . where rays travel from point A to point 5 . As r increases, the 
path from A io B becomes more haphazard due to the increased number 
of diffusers encountered. 

An analytical expression for the surface area of an ellipsoid exists, 
but contains a complete elliptic integral of the second kind, which has to be 
evaluated numerically. However, an approximate analytical form for the 
surface area can be derived, having a more intuitive form. The 
approximation for the surface area of the ellipsoid must hold in the limits 
t->cAB and'^r->oo. 

In the case of the surface area for the shortest reflection delay, 
which can equlvalently be viewed as the surface area for 45->oo, the 
reflectors lie on the surface of a cylinder of area: 

S = 27ibAB , (8) 

where b is the minor axis of the ellipsoid. 

Given that r is the reflected path length of a ray, then: 
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Therefore, for t -> cAB 0 , the number of reflectors is given by: 



N(r) = S = 7iAB^r^ 'AB^ , 



(10) 



In the case of the surface area for the longest reflection delay, 
which can equivalently be viewed as the- surface area for AB-^Q, the 
reflectors lie on the surface of a sphere, which has the simple form: 



where the path length r is twice the radius of the sphere. 
Hence, the approximate form for the number of diffusers as a 
function of time is given by: 



where c is the speed of light, t is the duration for the received 
reflected rays to travel from A io B and is the time for a ray to travel 
directly from A Xo B . 

So, over time the number of diffusers increases when the surface 
they lie on evolves from a cylinder into a sphere. Thus, the number of 
received rays in time, and therefore the received signal, also increase by 
the factor N{t) due to the Increased number of diffusers encountered by 
the rays. 

However, as time increases the amplitudes decay in free space, as 
mentioned in ill) above, by the amount: 



The effect of N{t) is reduced by free-space propagation decay to give the 
amount: 



(11) 




(12) 




(13) 




(14) 
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The effect of ttie three factors, i) to iii), of random-phase 
superposition, the numtier of diffusers, and free-space propagation decay, 
thereby combine to yield the following approximation for the diffuse 
componente of the received signal: 

(15) 



m = Ae 



AB 



+ r 



f(t) = Be^ 
where 



This can be simplified to: 



(16) 



5=^ and a^y^ 



The diffuse approximation can be viewed as a background function 
for the received signal r(0, where specular reflections represented* by the 
linear superposition signal basis can be superposed, resulting in the 
transformation of the MEDLL equation into the form: 



m=0 



+ 1 



(17) 



As the diffuse scattering process Is random, involves a large number of 
rays and exhibits an exponential decay, the above approximation provides 
a good characterisation of the diffuse process. 

The maximum likelihood estimator is given hereinbefore in Equation 

(1). 

The present invention hereby provides an improved estimator, with 
the capability to clean the correlation function of diffuse reflections, given 

by the following equation: 



.dx 
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(18) 



where s(t) is the replica signal, and B and a are free parameters. The 
other parameters all appear in the maximum likelihood estimator described 
hereinbefore. 



deriving from the approximation for the diffuse components of the received 
signal given in Equation (16), represents an improvement to the maximum 
likelihood (MEDLL) estimator. 

Figure 3 illustrates an example of the method by which the present 
invention may be implemented. 

It will be apparent to a person skilled in the art that this method 
could be embodied in either software or hardware, for example as part of 
the signal-processing unit of a radio receiver. 



The term 




(19) 
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CLAIMS 

1. Method of position determination in a radio system, the method 
comprising correlating (Step 1) a signal (R(t)) received at a unit (1) with a 
replica signal at the unit, and processing (Step 2) the correlated signal with 
an optimisation function comprising an exponential term in combination 
with a second term. ** " 

2. A method according to Claim 1 wherein the exponential term is 
in the form Be"°* (Step 2). 

3. A method according to Claim 1 or 2 wherein the second term is 
of the form: 

(Step 2). 

4. A method according to any preceding claim comprising effecting 
an Integration (Step 3) with the replica signal. 

5. A method according to any preceding claim comprising fitting the 
optimisation function and a Line-of Sight correlation function (Step. 4) with 
a set of parameters. 

6. A method according to Claim 5 comprising superposing the 
diffuse correlation output with a Line-of-Sight function output and fitting 
with correlation data of known values for the Line-of-Sight output. 

7. A method according to any preceding claim comprising first 
operating a multipath mitigation technique to effect correlation of the 
received and replica signals. 
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8. A method according to Claim 5 wherein the multipath mitigation 
technique comprises a Multipath Estimating Delay Locks Loop (MEDLL) 
technique (Step 1 ), 

9. A method according to Claim 5 wherein the multipath mitigation 
technique comprises a Minimum Mean Square Error (MMSE) technique. 

10. A computer program product directly loadable into the intemal 
memory of a digital computer, comprising software code portions for 
performing the method of any one or more of Claims 1 to 9 when said 
product is run on a computer. 

1 1 . A computer program directly loadable into the intemal memory 
of a digital computer, comprising software code portions for performing the 
method of any one or more of Claims 1 to 9 when said program is run on a 
computer. 

12. A carrier, which may comprise electronic signals, for a 
computer program of Claim 11. . 

1 3. Electronic distribution of a computer program product of Claim 
1 0 or a computer program of Claim 1 1 or a carrier of Claim 1 2, 

14. Apparatus for ppsltion determination of a radio system, the 
apparatus comprising means to correlate (13) a signal (R(t)) received at a 
unit (1) with a replica signal at the unit, and means (13) to process the 
correlated signal with a n optimisation function comprising an exponential 
term in combination with a second term. 

15. Apparatus according to Claim 14 wherein the exponential term 
is in the form Be'"*. 
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16. Apparatus according to Claim 14 or 15 wlierein tlie second 
term is of the form: 



- 17. Apparatus method according to any of Claims 14 to 16 
comprising means (15) to effect an integration with the replica signal. 

18. Apparatus according to any of Claims 14 to 17 comprising 
means (15) to fit the optimisation function and a Line-of Sight correlation 
function with a set of parameters. 

19. Apparatus according to Claim 18 comprising means (15) to 
superpose the diffuse correlation output with a Line-of-Sight function 
output and fit with correlation data of known values for the Line-of-Sight 
output. 

20. Apparatus according to any of Claims 14 to 19 comprising 
means to first operate a multipath mitigation technique to effect correlation 
of the received (R(t)) and replica signals. 

21. Apparatus -according to Claim 20 wherein the multipath 
mitigation technique comprises a IVIultipath E stlmating Delay Locks Loop 
(MEDLL) technique. 




22. Apparatus according to Claim 20 wherein the multipath 
mitigation technique comprises a Minimum Mean Square Error (MMSE) 
technique. * •"; 
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23. Method of position determination substantially as hereinbefore 
described w ith reference to, a nd/or a s illustrated I n, a ny o ne or more of 
the Figures of the accompanying drawings. 

24. Apparatus substantially as hereinbefore described with 
reference to, and/or as illustrated in, any one or more of the Figures of the 
accompanying drawings. 
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ABSTRACT 

METHOD OF POSITION DETERMINATION 

Unit 1 has antenna 10 to receive incoming signal R(t) and pass it to 
multiplier block 11 for multiplication with output C from the canier 
generator 12 in order to mix down signal R(t) which is passed to a block of 
correlators 13 for sequential correlation of the mixed down signal and the 
replica code from code generator, in order to produce correlation data 
outputs from the I (imaginary) and Q (real) channels using input signals 
from code generator 14. 

The correlation data is then p assed to signal processing block 1 5 , 
which adjusts appropriately both the code generator 13 and carrier 
generator 12 to provide accurate position measurement in indoor 
environments with multiple diffuse reflections. 



[Figure 2] 
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